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Abstract. We study the zero-point energy of a massless scalar field subject to 
s ! , spheroidal boundary conditions. Using the zeta-function method, the zero-point 

, energy is evaluated for small ellipticity. Axially symmetric vector fields are also 

ft \ • considered. The results are interpreted within the context of QCD flux tubes and 

p ' iln MIT ba.u, modi I. 
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1. Introduction 



It is well known that the zero-point energy of a quantum field is highly dependent 
on the geometry of any boundary conditions. Many geometries have been studied 
| including spherical, cylindrical and parallel plates £Q. These three geometries are 

the limiting cases of spheroidal geometry. By considering spheroidal geometry we 
may investigate mathematically how the zero-point energy changes as the geometry 
changes. In particular, the sign of the zero-point energy can be examined. 

QCD is often studied phenomenologically in terms of bag models [2]. The zero- 
■ point energy of a spherical bag has already been found Lattice QCD simulations 

have shown that there exist flux tubes between quark anti-quark pairs and also three 
^H. quark configurations 4 . We use spheroidal geometry to model these flux tubes. 

Using the infrared modified gluon bag model [5], it is found that the zero-point energy 
produces an attractive inter-quark force. 



' 2. Spheroidal geometry 

h ; 

The Cartesian equation for a spheroid is 

^ + 5 = i- ID 

a z b z 

If a < b (a > b), then the spheroid is prolate (oblate). If a = b, then the spheroid is 
a sphere of radius a. Prolate spheroidal curvilinear coordinates, (u,v,w), are related 
to Cartesian coordinates by 



x = ,/V(u 2 -l)(l-v 2 )cos(w) (2) 
y = fV(u 2 -l)(l-v 2 )sm(w) (3) 
z = fuv (4) 



where / = \/b 2 — a 2 and 1 < u < oo, — 1 < v < 1,0 < w < 2ir. In these coordinates, 
surfaces of constant u are confocal prolate spheroids. The prolate ellipticity, e, is a 
measure of deviation from a sphere and is defined by e = f/b £ [0, 1). 
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3. Scalar field 

3.1. Spheroidal junctions 

We consider a real, massless scalar field, <f>, confined within a prolate spheroid with 
Dirichlet boundary conditions. The field equation (with h = c = 1) and boundary 
condition are 

(V 2 +^ = 0, (5) 

0U=l/e=O, (6) 

where u) is the frequency. It can be shown that JSJ) is separable in prolate spheroidal 
coordinates [S]. Let <j) = U(u)V(v)W(w), then (J5J separates into three ordinary 
differential equations 

(5(<""- 1 >E)-* + A ' + 5^) t 'M-°- (7 » 

^((l-O+A-cV-^V-M.O, (8) 

^ + p')»'H = ». (9) 

where c = w/ and /i and A are separation constants. The solution to (|5J) is 
straightforward and requires fj, = m € Z. The solutions to J7J and JSJ that are 
regular and satisfy the boundary condition are the radial spheroidal function of the 
first kind, Rm]n(c, u), and the angular spheroidal function of the first kind, Sm]n(c, v), 
respectively. Here n is an integer such that n > \m\. In the spherical limit, 
reduces to the spherical Bessel differential equation and (JSJ reduces to the associated 
Legendre differential equation. In fact, the radial (angular) spheroidal functions can 
be expanded in a basis of spherical Bessel (associated Legendre) functions [Jj 

(u 2 - l)™/ 2 00 mn 

" 7 

K— — OO 

OO 

S { rkUc,v)= J2 <Q n (c)Pn+2k(v), (11) 
k— — oo 

where the coefficients a^ n {c) and d™^ n (c) satisfy certain three-term recurrence 
relations. Angular spheroidal functions are normalized in the same fashion as associate 
Legendre functions, whereas the radial spheroidal functions are normalized to have the 
same asymptotic form as spherical Bessel functions 

Rg] n {c,u) cu ^°° j n {cu). (12) 

The Dirichlet boundary condition can now be written in terms of radial spheroidal 
functions 

R£] n (z m , n , k e, 1/e) = 0, (13) 

where 

Zm,n,k — ^m,n,/c^ and ^ the k eigen- frequency. Note that b and e 

are used for the geometrical variables however, since only small ellipticity will be 
considered a and e could have equally been used. 
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3.2. Zero-point energy 

In terms of the eigen-frequcncies of a quantum field subject to some boundary 
condition, the zero-point energy is 



E 



7} ^2 U ra,n,k- (14) 



As expressed above, the zero-point energy is divergent and must be regularized. The 
prolate spheroidal zcta function is defined 



oo 



c m ,„( S )-M s ^^„, fc , (is) 

k=l 



where the z m ,n,fc come from <|13f) . M is an arbitrary scale with mass dimensions and the 
dependance on e has been suppressed. As written, l|15|) is only valid for Re s > 1. The 
regularized zero-point energy is defined in terms of a principal part prescription 8 , 



£7=i/xPP__ x 



^ ] Cm.n i s ) 



(16) 



Clearly, knowledge of the prolate spheroidal zeta function as s — > —1 is required. 
Fortunately, £(s) may be analytically continued to the region — 1 < Res < 0. Using 
the argument principal, l|15(l may be re-written 

Cm,n(*) = / dzz- 5 - 1 ln(i?W„(ze, 1/e)), (17) 

where the contour encloses all the zeros of Rm] n (ze, 1/e) . Because the radial spheroidal 
functions have the same asymptotic form as the spherical Bessel functions, we can 
follow an analogous argument to Romeo .9. , arriving at an integral expression for the 
spheroidal zeta function valid for —1 < Res < 0, 

Cm,«(«) = (&M) S ~ sin (^) j dxx-*- 1 In (2x exp{-x)T^ n (xe, 1/e)) , (18) 

where Tm]n(xe, 1/e) are the modified radial prolate spheroidal functions, modified in 
the same fashion as the modified spherical Bessel functions. If we make a Taylor 
expansion for small ellipticity, then 

£)&»,»(*) = (6/i) s -sin(y) X>^ s+1 / d.TX- s - 1 ln(2cxp(-^)s n (^)) 



m.n 7i—0 



n=0 17 V 



s„(ra) 



e 2 + O (e 4 ) , (19) 



where v — n + 1/2 and s n (i/x) are the Riccati-Bessel functions of imaginary argument. 
The first term in (|19fl is simply the spherical zeta function and has already been 
evaluated UJ; the second term is new and represents the first order correction to 
the spheroidal zeta function for small ellipticity. To evaluate the second term, the 
standard procedure is followed. Using the uniform asymptotic expansion of the 
modified spherical Bessel functions |l(Jj , the asymptotic form of the integrand is added 
and subtracted. The integrand is expanded to O (f~ 5 ), which is the same order as 
in The integral of the added terms can be expressed in terms of beta functions 
and the integral of the subtracted terms can be evaluated numerically. The sums over 
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n can be expressed in terms of Riemann zeta functions. The result, when Laurent 
expanded around s = — 1, is 

£ ^ = r» ( sL (ttt + ln( ^) + °- 0089 ) + (s + 1} 



1 / 2 



1 



_b[x \9457T \s + l 
Thus, using (fToj) . 



+ ln(&^) + 0.0023 + O (s + 1) 



e 2 + O (e 4 ) 



(20) 



26 V3157T 



ln(6/x) + 0.0089 



1 / 2 



■ ln(fyi) + 0.0023 



e 2 + O (e 4 



(21) 



2b V9457T 

where the superscript, I, signifies that this is the zero-point energy for internal eigen- 
frequencies only. 

If the field is confined to the outside of a prolate spheroid by Dirichlet boundary 
conditions, then the exterior eigen-frequencies are found. The exterior zero-point 
energy can be calculated using an analogous spheroidal zeta function. If internal and 
external eigen-frequencies are considered together, then the total zero-point energy is 

folate = ^ (0.0056 + 0.0019e 2 ) + O (e 4 ) . (22) 

As in the spherical case, there is no pole at s = — 1 in the total prolate spheroidal zeta 
function, so the principal part prescription is not needed. The total zero-point energy 
is independent of the mass scale. 

The zero-point energy of a real, massless scalar field subject to the same boundary 
conditions but on the surface of an oblate spheroid has also been calculated ^J. We 
quote here the result for the total energy 



E, 



oblate 



-(0.0056- 



0.0019e 2 ) + O (e 4 ) , 



where e 2 = a 2 /b 2 — 1 and is related to the oblate ellipticity. In both 
spherical part agrees with |12) . 



(23) 

and JUJ the 



3.3. Green's functions 

The results of the previous section can be confirmed using the Green's function 
approach. If all space is considered, then the zero-point energy is given by J3j 

E= lim i / — e-Ht-t') f dr 2u 2 g (r, r') , (24) 

where g (r,r') is the reduced Green's function, which satisfies 

(V 2 +uj 2 )g{r,r') = -S(r-r'). (25) 

The g (r, r') that satisfies the Dirichlet boundary condition on the surface of a prolate 
spheroid is 

g(r,r') = -iuj^2x mtn (c,v,w)X* ln (c,v',w') 



AR%{c,u)B^[ n {c,u') 



(i) 



<-m,n V 



u,u' < 1/e 
u,u' > 1/e 



(26) 
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Figure 1. Total zero-point energy of a scalar field subject to Dirichlct conditions 
on a prolate spheroidal (a) and an oblate spheroid (b) of small ellipticity. The 
dashed lines represent preliminary O (e 4 ) results. 



where the coefficients are 

A = B^ = ?^lIA (27 ) 

and Rm,n{c, u) are the radial prolate spheroidal functions of the third kind, which are 
given by equation IjlOjl with spherical Hankel functions of the first kind, hn\cu), as 
the basis functions (instead of j n (cu)). The prolate spheroidal harmonics are 

X m , n (c,v,w) = J \ +1 tl m l S £l( c > v ) eimW - ( 28 ) 
w 47r (n + my. 

If we expand l|24|) for small ellipticity and expand the resulting integrands 
asymptotically to O (v~ 2 ) the result agrees with to four decimal places. 

3.4- Results 

The results are shown in Figure ^ The preliminary O (e 4 ) results (dashed lines) 
reinforce that l|22ll and l|23[) are only valid for small ellipticity. 

It can be seen that 2bE pTO i atc increases as the sphere is deformed. This is expected 
as the limiting geometry is an infinite circular cylinder, for which it is known that the 
zero-point energy per unit length is positive 14 . 

In the oblate case, the zero-point energy becomes negative. Again, this is expected 
since the zero-point energy per unit area of two parallel plates is negative £Q. The 
precise point at which the zero-point energy is zero is outside the region where these 
results are valid. Work is currently being done to find this point. 

4. Vector field 

The vector Hclmholtz equation cannot be solved by separation of variables in prolate 
spheroidal coordinates [BJ. However, if the vector field is axially symmetric, then it 
can be separated into TE and TM modes. The TE and TM boundary conditions 



Zero-point energy in spheroidal geometries 



6 



for an axially symmetric field confined within a perfectly conducting prolate spheroid 
cavity are, respectively 



= l/e 



= 0, 



d_ 

du 



0. 



(29) 
(30) 



u—l/ e 



It should be noted that, due to the assumed axial symmetry, the radial spheroidal 
functions have m = I and now n G Z + . Adding the contributions from both the TE 
and TM modes and using l|16|) , the internal zero-point energy of an axially symmetric 
vector field is JT] 



16 



2b V 315tt 



ln(6/i) + 0.168 



— (0.037 ln(£^)- 0.057) 



O 



(31) 



In this calculation the integrands were expanded asymptotically to 0(y 2 ). The 
spherical part agrees with |To] . 



4.1. QCD 

Assuming the infrared modified gluon bag model 5] , the internal zero-point energy of 
a colour field is approximately 

^QCD ~ -8-Evj ( 32 ) 
where the factor of 8 comes from the eight gluon fields. The minus sign comes from 
modifying the free gluon propagator to one which vanishes in the infrared region. 

We model an inter-quark flux tube by a prolate spheroidal cavity with bag-like 
boundary conditions. Using the QCD zero-point energy is approximately 

£QCD«^(-l-34 + 0.45e 2 )+O(e 4 ), (33) 

where the scale independent part of l|31|l has been used and thus we must assume the 
colour fields are axially symmetric. If the flux tube is stretched, then (|33l) implies 
that the zero-point energy increases. Therefore, the zero-point energy contributes an 
attractive inter-quark force. The result also suggests that spherical bags are stable 
against deformations. 
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